How many of the scaling trends in pp collisions will be violated at -y/SNN = 14 TeV ? - 
Predictions from Monte Carlo quark-gluon string model. 
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Multiplicity, rapidity and transverse momentum distributions of hadrons produced both in inelas- 
tic and non-diffractive pp collisions at energies from y/s — 200 GeV to 14 TeV are studied within 
the Monte Carlo quark-gluon string model. Good agreement with the available experimental data 
is obtained and predictions are made for the collisions at top LHC energy y/s = 14 TeV. The model 
indicates that Feynman scaling and extended longitudinal scaling remain valid in the fragmentation 
regions, whereas strong violation of Feynman scaling is observed at midrapidity. The KNO scaling 
in multiplicity distributions is violated at LHC also. The origin of both maintenance and violation 
of the scaling trends is traced to short range correlations of particles in the strings and interplay 
between the multi-string processes at ultra-relativistic energies. 
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I. INTRODUCTION 



The recent interest to general features of elementary 
hadronic interactions, especially to characteristics of pp 
collisions, at ultrarelativistic energies is manifold. First 
of all, these collisions are conventionally used as refer- 
ence ones to reveal the nuclear matter effects, such as 
strangeness enhancement, nuclear shadowing, collective 
flow etc., attributed to formation of a pattern of hot 
and dense nuclear matter and the quark-gluon plasma 
(QGP) in the course of ultrarelativistic heavy-ion colli- 
sions (see [l| and references therein). Although the for- 
mation of the QGP and/or collective behavior was not 
found yet in pp collisions at energies up to the Tevatron 
energy y/s — 1.8 TeV, physicists are discussing the pos- 
sibility to observe, e.g., elliptic flow in pp interactions 
at y/s = 7 TeV accessible for the Large Hadron Col- 
lider (LHC) at CERN at present. This limit may be 
raised to y/s = 14 TeV in the nearest future. Because of 
the huge amount of energy deposited in the overlapping 
region the pp systems might be similar to A + A colli- 
sions at nonzero impact parameter at lower energies @, Q 
and, therefore, demonstrate collectivity. An alternative 
approach developed in [|| considers the flow effects in 
hadronic interactions as initial state effects linked to cor- 
relation between the transverse momentum and position 
in the transverse plane of a parton in a hadron. This 
problem should definitely be clarified in future. Then, 
the problem of multiparticle production in elementary 



hadronic collisions is not fully solved yet. Here for hard 
processes with large momentum transfer the running cou- 
pling constant as is small, that allows for application 
of the perturbative quantum chromodynamics (QCD). 
For soft processes with small momentum transfer, which 
give dominant contribution to high energy hadronic in- 
teractions, the as is close to unity and, therefore, non- 
perturbative methods should be applied. Many micro- 
scopic models I5l 4l3l based on the string picture of par- 
ticle production [14[ have been successfully employed to 
describe gross features of hadronic collisions at relativis- 
tic and ultrarelativistic energies, whereas the statistical 
approach pioneered more than 50 years ago by Fermi and 
Landau [ll, [f6| is not ruled out. To make predictions for 
the LHC in the latter case one has to extrapolate the 
data obtained at lower energies to the high energy re- 
gion. It was found quite long ago that, despite the com- 
plexity of a reaction with tens or more particles in a final 
state, multiparticle production in pp collisions exhibits 
several universal trends, such as (In y/s) 2 dependence of 
total charged particle multiplicity [13], Feynman scaling 
[l8j and related to it extended longitudinal scaling 
KNO scaling [2(| and so on. Similar trends were found 
later on in proton-nucleus and nucleus-nucleus collisions 
as well (for review see, e.g., [HI). On the other hand, 
the description of ultrarelativistic hadronic interactions 
in the framework of Color Glass Condensate (CGC) the- 
ory [iH leads to a universal power-law behavior of, e.g., 
density of produced charged particles per unit of rapidity 
and their transverse momentum 23 1 . 
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The aim of the present article is to study the main 
characteristics of pp interactions at energies from y/s = 
200 GeV to top LHC energy y/s = 14 TeV. We employ 
the Monte Carlo (MC) realization [7j of the quark-gluon 
string model (QGSM) [24| based on Gribov's Reggeon 
field theory (RFT) j25( that obeys both analyticity and 
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unitarity requirements. The features of the model are 
described in Sec. [TT] in details. Comparisons with avail- 
able experimental data for pp and pp collisions at energies 
y/s > 200 GeV, including the recently measured pp inter- 
actions at \fs = 2.36 TeV and \fs — 7TeV, and predic- 
tions for the top LHC energy yfs — 14 TeV are presented 
in Sec. Mil Here exclusive contributions of soft and hard 
processes to particle rapidity and transverse momentum 
spectra are studied. Special attention is given to the ori- 
gin of violation of the KNO-scaling, violation of the Feyn- 
man scaling at midrapidity and its maintenance in the 
fragmentation region. Obtained QGSM results are also 
confronted to the predictions of other microscopic and 
macroscopic models. Conclusions are drawn in Sec. IIVI 



II. 



QUARK-GLUON STRING MODEL AND ITS 
MONTE CARLO REALIZATION 
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As was mentioned in the Introduction, the description 
of soft hadronic processes cannot be done within the per- 
turbativc QCD. Therefore, the quark-gluon strin g m odel 
24] employs the so-called 1/N series expansion |26l. I27I 
of the amplitude for processes in QCD, where N is either 
number of colors N c [2(| or number of flavors Nf [27} ■ It 
this approach the amplitude of a hadronic process is rep- 
resented as a sum over diagrams of various topologies, so 
the method is often called topological expansion. It ap- 
pears that at high energies and small momentum trans- 
fer the arising diagrams are similar [28|, [29[ to processes 
describing the exchange of Regge singularities in the t- 
channel. For instance, planar diagrams correspond to the 
exchange of Reggeons, and cylinder diagrams correspond 
to reactions without quantum number exchange in the 
i-channel, i.e., taking place via the Pomeron exchange, 
where Pomeron is a composite state of the reggeized glu- 
ons. Processes with many- Pomeron or many-Reggeon ex- 
changes are also possible. To find the amplitude of mul- 
tiparticle production one has to cut the diagrams in the 
s-channel, and the physical picture of quark-gluon strings 
arises. Namely, new particles are produced through the 
formation and break-up of quark-gluon strings or exited 
objects consisting of quarks, diquarks and their antistates 
connected by a gluon string. 

Figure [1] shows the subprocesses with particle creation 
taken into account in the current Monte Carlo version of 
the QGSM [7| for pp collisions at ultrarelativistic ener- 
gies. The inelastic cross section consists of three terms 

a inel( S ) = a P( S ) + °"Sd(s) + <JDd(s) , (1) 

where erp(s) is the cross section for the multi-chain pro- 
cesses described by the cylinder diagram and diagrams 
with multi- Pomeron scattering [Fig.QTa)], &sd(s) by the 
single diffractive processes with small [Fig. [ljb)] and 
large [Fig. QJc)] mass excitation, corresponding to the 
triple- Reggeon and triple- Pomeron limit, respectively, 
and (Tfljj(s) by the double diffractive diagram [Fig.QJd)]. 



FIG. 1: Diagrams of particle production processes included 
in the modeling of pp interactions at ultrarelativistic energies. 
See text for details. 

Other diagrams that are relevant at low and intermedi- 
ate energies, such as undeveloped cylinder diagram or 
diagram with quark rearrangement [7| , play a minor role 
here because their cross sections rapidly drop with rising 
s. The statistical weight of each subprocess is expressed 
in terms of the interaction cross section for the given sub- 
process CTj(s) 

Wj = cn(s)/a ine i{s) . (2) 

Then, the hadron inelastic interaction cross section 
o~inei{s) — <7tot{s) — cr e i(s) is split into the cross section 
for single diffractive interactions <jsd(s) and the cross 
section for non-diffractive reactions o~nd(s), as it is usu- 
ally done in analysis of experimental data. By means of 
the Abramovskii-Gribov-Kancheli (AGK) cutting rules 
[30| the inelastic non-diffractive interaction cross section 
o~n d (s) can be expressed via the sum of the cross sections 
for the production of n = 1,2,... pairs of quark-gluon 
strings, or cut Pomerons, and the cross section of double 
diffractive process 

oo 

o-nd(s) = <r n (s) + (J DD (s) . (3) 

n=l 

To find <7„(s) one can rely on the quasi-eikonal model 
[2^ . l33j which states that 

00 

(T to t(s) = ^2<?n(s) = (Tpf (|) , (4) 
I n— 1 & \ 

an{s) = ^( 1 " exp( ~ 2,) S^r) * ^ 
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The cross section do corresponds to diffraction contribu- 
tion. The parameters 7p and Rp are Pomeron-nucieon 
vertex parameters, quantity A = ap(0) — 1, and ap(0) 
and a'p is the intercept and the siope of the Pomeron 
trajectory, respectively. The quantity C takes into ac- 
count the deviation from the pure eikonal approximation 
(C = 1) due to intermediate inelastic diffr active states, 
£ = In (s/sq) and so is a scale parameter. 

For the diffractive processes displayed in Fig. [IJb) and 
Fig. [TJc) the fractions of momenta of initial hadrons car- 
ried by the sea quark pairs x q q are determined according 
to distribution 



U qq( X qq) « 



(1 



< (/(/ 



with the cut for the minimum values 

0.14 
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(10) 
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The transverse momentum distribution of (anti)quarks 
in a proton in the low-mass excitation process shown in 
Fig. DIb) is given by 



fq(PT)dp T 



exp {-hip?) dpr , 



(12) 



where the slope parameter b\ = 20(GeV/c)~ 2 . Then, 
it is assumed that the valence (anti)diquark in the 
(anti)proton carries a transverse momentum equal in 
magnitude and opposite in sign to the sum of trans- 
verse momenta of the other (anti)quarks. The number of 
quark-gluon strings increases with collision energy, thus 
the average transverse momentum of the (anti)diquark 
rises also. 

The quantitative description of single diffractive and 
double diffractive processes at high energies was done 
in QGSM recently [HI, HI] in terms of dressed tripplc- 
Reggeon and loop diagrams. The results obtained in [31j 
for the cross sections of the diffractive processes are uti- 
lized in our MC model via the parametrizations 

<Jsd{s) = 0.68 (l + y) In (0.6 + 0.2 s) , (13) 



o-dd(s) 



1.65 + 0.27 Ins 



(14) 



Soft processes dominate the particle production in 
hadronic interactions at intermediate energies. With 



the rise of the collision energy hard processes result- 
ing to formation of hadronic jets with large transverse 
momenta pt become important also. To take into ac- 
count the jet formation and, on the other hand, to de- 
scribe simultaneously the increase of the total and in- 
elastic hadronic interaction cross section with rising y/s 
the eikonal model was properly modified in [34j by intro- 
ducing the new term that represents the hard Pomeron 
exchange. The cut of the hard Pomeron leads to for- 
mation of two hadronic jets, see Fig. [5J Therefore, the 
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FIG. 2: String formation in hard gluon-gluon scattering and 
soft Pomeron exchange in proton-proton collision. 

eikonal u(s, b), that depends on the center-of-mass en- 
ergy y/s and the impact parameter b, can be decomposed 
onto the terms corresponding to soft and hard Pomeron 
exchange: 



u(s, b) = u so f t (s, b) + U h ardi.S, b) . 



(15) 



The inelastic hadronic cross section Oi ne \ (s) is connected 
to the real part of the eikonal u R (s, b) as 



OO 

<Tinel(s) = 2tt J {1 - exp [-2u R (s : b)] } bdb 



(16) 



Recall that the concept of (semi)hard Pomeron is nowa- 
days a common feature of all RFT-based MC mod- 
els [1, [H, [T3, [H, [I!] designed for the description of 
hadronic and nuclear interactions at ultrarelativistic en- 
ergies. Other microscopic MC models also rely on the 
picture of soft+hard eikonal approach [10( . 

Following [3J, [35| , both soft and hard eikonal can be 
expressed as 



U foft/hard( S i b ) = Z S oft/hard(s) exp 

where [cf. Eqs. Q-©] 
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^^soft/hard( s ) 



z soft/hard{s) 



IP 



^soft/hard{ s ) \ s 



a P (0)-l 



(17) 



(18) 



Koft/hard{s) = R P + a' P In ( — ] . (19) 
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TABLE I: Parameters of the soft and hard Pomeron used in 
the current version of the QGSM. 



Parameter 


Soft Pomeron 


Hard Pomeron 


a P (0) 
a' P 
~fp 
Rp 


1.1855 
0.175 
0.832 
2.6 


1.3217 


0.021 
2.4 



Numerical values of the slopes and intercepts of the 
Pomeron trajectories and parameters of the hadron cou- 
pling to the Pomeron used in the model fit to experi- 
mental data are listed in Table Q] Note that these values 
deviate from the parameters of the soft and the hard 
Pomeron obtained in [34|, [53] from the cross section of 
minijets measured by the UA1 Collaboration. To de- 
scribe the LHC data at energies above y/s = 900 GeV it 
was necessary to increase the soft Pomeron intercept to 
ap(0) — 1 « 0.185 and simultaneously to reduce the value 
of its slope parameter a' P to 0.175. 

Then, the AGK cutting rules enable one to express the 
inelastic cross section as 



where 



(Jinel(s) = ^ U V ( S ) > 
i,j—0;i-\-j>l 



i(s) = 2n J bdb exp [-2u R {s, b)] 
o 

^ 2u« /t ( S ,6)]' [2u« rd (s,&)] : 



(20) 



(21) 



The last expression can be used to determine the num- 
ber of quark-gluon strings and hard jets via the number 
of cut soft and hard Pomerons, respectively. At very 
high energies one has to take into account the effects of 
shadowing of partonic distributions both in nucleons and 
in nuclei. In the Reggeon calculus such processes corre- 
spond to the so-called enhanced diagrams 38] describing 
the interactions between Pomerons. These diagrams are 
not implemented yet in the current MC version of the 
QGSM. 

It is worth noting here that the multi-Pomeron ex- 
changes become very important with increasing cm. en- 
ergy of hadronic collision. For instance, the contribution 
of a single-cylinder diagram to the scattering amplitude 
is proportional to (s/s ) apl - Q ^~ 1 , ap(0) > 0. In contrast, 
the contributions coming from the n-Pomeron exchanges 
grow as (s/so) nA . Although in the framework of the 
1/A^-expansion the n- Pomeron exchange amplitudes are 
suppressed by factor 1/N 2n , the quickly rising term s nA 
dominates over the suppression factor at ultrarelativistic 
energies. 

There is no unique theoretical prescription for mod- 
eling the fragmentation of a string with a given mass, 



TABLE II: Total, elastic, multi-chain, single diffraction and 
double diffraction cross sections of pp collisions calculated by 
the QGSM. 

y/s (GeV) a tot (mb) a e \ (mb) a P (mb) o S d (mb) a DD (mb) 



200 


51.62 


9.67 


31.12 


6.12 


4.51 


546 


60.83 


12.51 


35.72 


7.48 


5.05 


630 


62.25 


12.97 


36.42 


7.67 


5.13 


900 


65.85 


14.15 


38.19 


8.16 


5.32 


1800 


72.97 


16.55 


41.61 


9.10 


5.70 


2360 


75.74 


17.50 


42.92 


9.47 


5.84 


7000 


86.60 


21.31 


47.91 


10.95 


6.43 


14000 


93.07 


23.61 


50.76 


11.89 


6.80 



momentum and quark content into hadrons. In the 
presented version of the QGSM the Field-Feynman al- 
gorithm [3!| is employed. It enables one to consider 
emission of hadrons from both ends of the string with 
equal probabilities. The break-up procedure invokes the 
energy-momentum conservation and the preservation of 
the quark numbers. The transverse momentum of the 
(di) quarks in the vacuum pair is determined by the 
power-law probability 



f(pr) dvr 



3Db 2 (s) 
7r(l + Df4) 4 
0.325 + 0.016 Ins 



(22) 
(23) 



with D = 0.34 (GeV/c)- 2 . 

Further details of the MC version of QGSM and its 
extension to h + A and A + A collisions can be found in 

ana Us- 



ui. COMPARISON WITH DATA AND 
PREDICTIONS FOR LHC 

For the comparison with model results concerning the 
pseudorapidity and transverse momentum distributions 
we used experimental data obtained by the UA5 Collab- 
oration for proton - antiproton collisions at cm. energies 
y/s = 200 GeV, 546 GeV and 900 GeV [H, by the UA1 
Collaboration for pp collisions at y/s = 546 GeV [43| , by 
the CDF and the E735 Collaborations for pp collisions 
at Vs = 1800 GeV 0,113, and recent CERN LHC data 
obtained in pp interactions at y/s = 900 GeV. 2360 GeV, 
and 7TeV by the ALICE Collaboration gjg| and by 
the CMS Collaboration [5(| ■ At such high energies 
the annihilation cross section is almost zero and the main 
characteristics of particle production in pp and pp inter- 
actions are essentially similar. 

Total and elastic cross sections are listed in Table [H] 
together with the cross sections of multi-chain, single- 
and double-diffraction processes for energies ranging from 
y/s = 200 GeV to y/s = 14TcV. Compared to those 
at y/s = 900 GeV, a tot , a e i and osd increase at y/s = 
14TeV by nearly 50%, whereas <tdd increases by less 
than 30%. 
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□ UA5 200 GeV O CDF 1.8 TeV 
❖ UA1 546 GeV * CMS 2.36 TeV 
A UA5 900 GeV * CMS 7 TeV 
lines - QGSM 




4 6 

p T (GeV/c) 

FIG. 3: (Color online) Transverse momentum distributions 
of the invariant cross section of charged particles in NSD pp 
collisions obtained in QGSM at \y\ < 2.5 for all energies in 
question. Experimental data are taken from [42T4441 . [501] . 



The transverse momentum distributions of the invari- 
d 3 a 

ant cross section E—^ divided to a tot for charged parti- 
al 

cles in non-single diffractive (NSD) pp collisions at all en- 
ergies in question are presented in Fig. [3l We see that the 
QGSM reproduces the experimental data in a broad en- 
ergy range pretty well. The spectra become harder with 
increasing y/s, thus the average transverse momentum 
of produced hadrons should grow also. Fig. |4] displays 
the (pt) of charged particles in NSD pp events calcu- 
lated in QGSM and compared to experimental data. We 
assume here 5% systematic errors for the extraction of 
mean pt because we do not apply any extrapolation pro- 
cedure to the generated spectra, as it is usually done in 
the experiments. Results of the fit of model simulations 
to quadratic logarithmic dependence and to power-law 
dependence are as follows 

(p T ) = 0.42 - 0.006 In s + 0.0022 In 2 s , 
(p T ) = 0.27 + 0.212£ 0115 . 

In the last expression E = y/s/ 2 and the exponent 0.115 
is not a free parameter. According to (23j, this exponent 
is just half of the exponent of power-law fit to dN/di] 
distribution (see below). As one can see in Fig. |H the 
difference between two paramctrizations of mean pt is 
negligible even for top LHC energy y/s = 14 TeV. 

To study the interplay between the soft and hard 
processes we show separately in Fig. [5] their fractional 
contributions and combined results for pp collisions at 
y/s = 900 GeV, 2.36 TeV, 7 TeV and at top LHC energy 
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FIG. 4: (Color online) Average transverse momentum as 
function of y/s. Squares present the QGSM results, other 
symbols denote experimental data from 42-45, 48, 50]. Solid 
and dashed lines are fit to logarithmic and power-law depen- 
dences, respectively. See text for details. 



y/s = 14 TeV. Moreover, the pT-dependence of the un- 
derlying soft processes from the collisions with at least 
one hard Pomeron is displayed in these plots as well to- 
gether with the experimental data of the CMS Collabo- 
ration. As expected, the soft processes dominate at low 
and intermediate transverse momenta, whereas at higher 
transverse momenta the major contribution to the cross 
section comes from the minijets. The crossover between 
the hard and soft branches takes place at pt ~ 2.8 GeV/c 
for the reactions at y/s = 900 GeV. It is shifted to 
p T w 2.2GeV/c at y/1 = 14 TeV. The slopes of the p T - 
spectra for both soft and underlying soft processes are 
similar. At y/s — 7 TeV and 14 TeV both lines coincide, 
i.e., the contributions to the invariant cross sections from 
barely soft Pomeron processes are equal to those from the 
soft Pomerons exchanges, accompanied by one or more 
hard Pomeron ones. 

Let us briefly recall the main assumptions and pre- 
dictions of the hypothesis of Feynman scaling [l8j]. It 
requires scaling behavior of particle spectra within the 
whole kincmatically allowed region of the Feynman scal- 
ing variable xf = p\\/p™ ax or, alternatively, cm. rapidity 
y* at ultrarelativistic energies s —¥ oo. In addition, the 
existence of non- vanishing central area |xf| < %o , ~ 
0.1 is postulated. In terms of rapidity this central region 
increases with rising y/s as 



(Ay* 



2 In [ajo 



s/m T 



(24) 



provided the transverse mass vtit — \J ttlq + p T is finite. 
In contrast, the fragmentation region remains constant 



{Ay*) f rag « ln(l/a; Q ) 



(25) 



From here it follows that (i) in the central area the par- 
ticle density p C ent(y* ,Pt, s) depends on neither y* nor 
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FIG. 6: (Color online) The charged particle pseudorapidity 
spectra for (a) inelastic and (b) non-single-diffractive events 
calculated in QGSM for pp collisions at \fs = 200 GeV (6), 
546 GeV (5), 900 GeV (4), 2.36TeV_(3), 7TeV (2) and 14 TeV 
(1). Data are taken from S3, S3, ■ 



FIG. 5: (Color online) Transverse momentum distributions 
of the invariant cross section of charged particles in NSD pp 
collisions at (a) yfs = 900 GeV, (b) yfs = 2360 GeV, (c) yfs = 
7 TeV and (d) = 14 TeV calculated in QGSM. Combined 
contribution of all processes and, separately, of only soft, hard 
and underlying soft sub-processes are shown by solid, dash- 
dotted, dashed and dotted lines, respectively (see text for 
details). Experimental data plotted in (a), (b) and (c) panels 
are taken from W9, [501 . 



y/s, i.e. p cen t = Pcent{pr), and rapidity spectra of parti- 
cles have, therefore, a broad plateau; (ii) this area gives 
a main contribution to average multiplicity of produced 
hadrons; (iii) contribution to the average multiplicity 
from the fragmentation regions is energy independent. 



The charged particle pseudorapidity 

1 dcf ine i 1 da NSD 



and 



spectra 
for inelastic and 



Cinei drj <Jnsd drj 

non-single-diffractive events, respectively, are displayed 
in Fig. EJa) and Fig. [U(b) together with the pp(pp) 
data at ^ = 200 GeV, 546 GeV, 900 GeV, 2.36 TeV 
and 7 TeV. QGSM predictions for yjl = 14 TeV are 
plotted here also. The model gives good description of 
these distributions within the indicated energy range 
except, maybe, not very distinct dip at midrapidity for 
the lowest energy in question \fs = 200 GeV. For pp 
collisions at top LHC energy QGSM predicts further 



increase of the central particle densities to 



dNi, 



dy 



|t?=o 



6.2 



dN 



NSD 



I r,=0 



7.1 



Compared to the y/s = 7 TeV, the rise of the central 
particle density at y/s = 14 TeV is expected to be about 
25%. 

In Fig. [7] the charged particle density at rj = is pre- 
sented as a function of the cm. energy y/s for inelas- 
tic (upper plot) and non-single diffractive (bottom plot) 
events. The experimental data for inelastic collisions be- 
low yfs — 546 GeV are well described by a linear depen- 
dence on Ins [42]. The striking evidence of first LHC 
results for pp interactions at \/s = 900 GeV, 2.36 GeV 
and 7 TeV is the quadratic dependence of the increase 
of midrapidity density of charged particles with rising 
Ins The theory of CGC suggests power- law rise 

[23j . In the QGSM these trends hold also, and the fit- 
ting parametrizations for cm. energies from 200 GeV to 
14 TeV are 



dN, 



inel 



drj 



dN 



NSD 



di] 



=0 (s) = 5.87- 0.74 Ins + 0.39 In 2 s 



v =o(s) = 6. 58 - 0.84 Ins + 0.45 In 2 s , 



dN NS p 
drj 



„=o(s) = 0.77£ 



0.23 
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FIG. 7: (Color online) The charged particle density at midra- 
pidity as a function of yfs for (a) inelastic and (b) non-single 
diffractive collisions. Dashed lines show the results of the 
fit to expression a + b Ins + c In 2 s, dotted line (b) - to the 
power- law dependence d y/s . 



FIG. 8: (Color online) The distributions 



do N SD 



as 



unsd dy 

functions of rapidity difference y—y ma x obtained in QGSM for 
energies yfs = 200 GeV, 546 GeV, 900 GeV, 1.8 TeV, 2.36 TeV, 
7 TeV and 14 TeV. 



As in the mean-pr case, there is a hairwidth difference 
between the two curves representing the logarithmic and 
the power-law fit, respectively. Indicating further in- 
crease of particle density at r\ — with rising energy, the 
model favors violation of the Feynman scaling at midra- 
pidity, otherwise the particle density there should not 
depend on y/s. 

It is interesting to compare the QGSM predictions for 
the charged particle multiplicity in pp collisions at LHC 
with that obtained by the extrapolation of pseudora- 
pidity distributions measured at lower energies. This 
method [2l| employs the energy independence of the 
slopes of pscudorapidity spectra combined with logarith- 
mic proportionality to ^fs of both the width and the 
height of the distributions. Therefore, any data set from 
Fig.JSJa) and Fig.ETb) can be used for the extrapolation, 
and the results are 12111 



dN 



NSD 



dl) 



\v=o 



4.6 ±0.4 



dN 



NSD 



dl] 



v=±2 = 5.25 ±0.7 



These predictions are significantly lower than the recent 
experimental data from LHC and the QGSM calcula- 
tions. 

Another feature that is closely related to Feynman scal- 
ing is the so-called extended longitudinal scaling [l9[ ex- 
hibited by the slopes of (pseudo)rapidity spectra. In the 
QGSM these slopes arc identical in the fragmentation re- 
gion ybeam > —2.5 as shown in Fig. [HI where the distribu- 



tions 



1 da N sD 



ctnsd 



dy 



are expressed as functions y — y n 



QGSM indicates that the extended longitudinal scaling 
remains certainly valid at LHC. This result contradicts 
to the prediction based on the statistical thermal model 
[5l| . The latter fits the measured rapidity distributions 
to the Gaussian, extracts the widths of the Gaussians 
and implements the energy dependence of the obtained 
widths to simulate the rapidity spectra at LHC. The ex- 
trapolated distribution was found to be much narrower 
51] compared to that presented in Fig. [5J We are ea- 
gerly awaiting the LHC measurements of pp collisions in 
the fragmentation regions to resolve the obvious discrep- 
ancy. Note, that experimentally the extended longitudi- 
nal scaling was found to hold to 10% in a broad energy 
range from ^/l = 30.8 GeV to 900 GeV [13]. 

The emergence of the extended longitudinal scaling as 
well as Feynman scaling in the QGSM is not accidental. 
It arises due to short range correlations in rapidity space. 
The correlation function of particle i and particle j, pro- 
duced in the string fragmentation, drops exponentially 
with rising rapidity difference 



C{yi,Vj) 



d 2 a 



da 



da 



Vinel dyidyj a ine i dyi a me i dyj 
cx exp[-X(yi-yj)} , (26) 



and, therefore, the particles with large rapidity difference 
are uncorrelated. Consider now the inclusive process 1 + 
2 — > i + X. Its single particle inclusive cross section 



fi = E 



d 3 a t _ d 2 a{yi - y t , y l - y2,pj T ) 



d 3 p 



dyid 2 pi 



(27) 
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becomes independent of yi — 2/2 at sufficiently high col- 
lision energy in the fragmentation region of particle 1, 
provided y\ - yi « 1 and yi - y 2 « y\ - yi > 1. Thus, 
the inclusive densities rii = fi/umel are determined by 
only two variables 

Hi = <p(yi - yi, pj T ) . (28) 
Recalling that the Feynman variable 

x i f = « exp [- (yi - w)] , (29) 
p \\ 

one arrives from Eq. (l23[) to the condition of Feynman 
scaling 

rn=il;{xf,p t iT ). (30) 

The invariant distribution F(xf) = 

— F f Ecm , ~ t 9 rfpr is displayed in Fig. |9] for 

all charged particles from the collisions at energies 
from = 200 GeV, to 14TeV. The scaling seems to 
hold within 20% of accuracy in the fragmentation region 
at 0.1 < x F < 0.2 only. 



energy independent function. KNO-scaling was found to 
hold up to ISR energies, y/s < 62 GeV, despite the ap- 
parent failure of the Feynman scaling hypothesis in the 
central region |xf| < xo- Violation of KNO-scaling was 
predicted within the QGSM in [24j. Later on the viola- 
tion was observed experimentally by the UA5 and UA1 
collaborations in pp collisions at y/s = 546 GeV [42j . The 
origin of this phenomenon in the model is the follow- 
ing. At ultrarelativistic energies the main contribution 
to particle multiplicity comes from the cut-Pomerons, 
and each cut results to formation of two strings. Short 
range correlations inside a single string lead to a Poisson- 
like multiplicity distribution of produced secondaries. At 
energies below 100 GeV the multi-string (or chain) pro- 
cesses are not very abundant and invariant masses of the 
strings are not very large. Therefore, different contribu- 
tions to particle multiplicity overlap strongly, and KNO- 
scaling is nearly fulfilled. With rising yfs the number of 
strings increases as (s/so) A and their invariant masses 
increase as well. This leads to enhancement of high mul- 
tiplicities, deviation of the multiplicity distribution from 
the Poisson-like behavior and violation of KNO-scaling 

am. 
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FIG. 9: (Color online) The invariant distribution F(xf) in 
non single diffractive pp collisions obtained in QGSM at yfs = 
200 GeV, 546 GeV, 900 GeV, 1.8 TeV, 2.36 TeV and 14TeV. 

Another scaling dependence is known as Koba-Nielsen- 
Olesen or KNO-scaling [2(|. Initially it was also derived 
from the hypothesis of Feynman scaling, but later on it 
appeared that both hypotheses are of independent origin. 
The KNO-scaling claims that at yfs — > 00 the normalized 
multiplicity distribution just scales up as In s or, equiva- 
lently, that 



(n) cr n 



= * 



(31) 



with a n being the partial cross section for n-particle pro- 
duction, (n) - the average multiplicity and ^(n/(n)) - 
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ALICE pp @ 2360 GeV lr]l < 1.3 *400 - 




FIG. 10: (Color online) Charged-particle multiplicity distri- 
butions in \rj\ < 0.5, \r)\ < 1.0 and \rj\ < 1.3 intervals, obtained 
in QGSM for pp collisions at y/s = 900 GeV (upper plot) and 
at yfs = 2360 GeV (bottom plot). Open symbols show the 
corresponding ALICE data [471 ]. 
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Before studying the violation of KNO-scaling at LHC, 
we compare in Fig. QJJ] the QGSM calculations with the 
ALICE data. In this figure the multiplicity distribu- 
tions of charged particles calculated in NSD pp events 
at yfs = 900 GeV and </s = 2.36 TeV in three central 
pseudorapidity intervals are plotted onto the experimen- 
tal data. The agreement between the model results and 
the data is good. Moreover, the QGSM demonstrates a 
kind of a wavy structure mentioned in [47j • As we will see 
below, such a wavy behavior in the model can be linked 
to nrocesses ening- via. the manv-Pomeron exchanges. 
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FIG. 11: (Color online) Charged-particle multiplicity dis- 
tributions in the KNO-variables in QGSM nondiffractive pp 
collisions at %/s = 546 GeV, 900 GeV, 1.8 TeV and 14 TeV. 
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FIG. 12: (Color online) Charged-particle multiplicity dis- 
tribution (solid line) for processes going via the exchange 
of 7i > 1 soft Pomerons in pp collisions at \/s = 14 TeV. 
Contributions of first four terms are shown by dash-dotted 
(n — 1), double-dash-dotted (n = 2), dashed (n — 3) and 
dotted (n = 4) lines, respectively. 



The multiplicity distributions of charged particles ob- 
tained in QGSM for NSD pp collisions at all ener- 
gies in question is presented in Fig. [TTJ for the interval 
|?7| < 2.4. Although the differences between the neigh- 
bor energies seem not to be very dramatic, the ten- 
dency in the modification of the distributions is quite 
clear. - The high-multiplicity tail is pushed up, maxi- 
mum of the distribution is shifted towards small values of 
n c h/(n c h) and the characteristic "shoulder" in the spec- 
trum becomes quite distinct, as presented by the distri- 
bution for top LHC energy. Another interesting obser- 
vation is the unique intersection point for all distribu- 
tions. All curves cross each other at z ~ 2.3 as can be 
clearly seen in the bottom plot of Fig. [Til where the ra- 
tio (n ch )P(z)\ 54:6GeV /(n ch )P(z)\ 7TeV is displayed. Note 
that the aforementioned pseudorapidity range \i]\ < 2.4 
is not sufficient to observe the multi-hump structure in 
the KNO-plot predicted in [29j (see also [Hj]) for the 
full phase-space. To clarify the role of multi-Pomeron 
processes in violation of KNO-scaling explicitly, Fig. fl2l 
shows the contribution to particle multiplicity diagram 
coming from the processes with different number of soft 
Pomerons in pp collisions at y/s = 14 TeV. The maxima 
of distributions for multi-Pomeron processes are moved 
in the direction of high multiplicities thus lifting the 
high-multiplicity tail. The pronounced peak in the low- 



multiplicity interval arises solely due to single Pomeron 
exchange. 

IV. CONCLUSIONS 

Multiplicity, transverse momentum and 
(pseudo)rapidity distributions of hadrons produced 
in pp interactions at energies from y/s — 200 GeV to 
14 TeV are studied within the Monte Carlo quark-gluon 
string model. Parameters of soft and hard Pomerons are 
determined from the fit to recent LHC pp data. Com- 
pared to the fit to lower energies, it was found necessary 
both to increase the intercept of soft Pomeron and to 
reduce its slope parameter. Other parameters, such as 
total cross sections, cross sections of single diffr active 
and double diffractive processes etc, are taken from 
theoretical considerations. The model simulations of 
pseudorapidity, transverse momentum and multiplicity 
spectra of secondaries are in a good agreement with 
the corresponding experimental data obtained in pp 
and pp collisions at Tevatron and at CERN energies. 
Predictions are made for pp interactions at top LHC 
energy y/s — 14 TeV. It is shown that within the 
examined energy range one cannot distinguish between 
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the standard" logarithmic dependence (oc In 2 s) and 
novel power- law approximation (oc E x ), employed for 
particle densities and for their mean px in Regge theory 
and in theory of Color Glass Condensate, respectively. 

Several scaling properties observed in particle produc- 
tion at relativistic energies have been examined. QGSM 
favors violation of Feynman scaling in the central rapid- 
ity region and its preservation in the fragmentation ar- 
eas. Extended longitudinal scaling is shown to hold at 
LHC. This scaling is also attributed to heavy-ion colli- 
sions at energies up to ^/s = 200 GeV. Extrapolations 
based on statistical thermal model predict its violation 
at LHC, thus implying vanishing of Feynman scaling for 
nuclear collisions in the fragmentation regions as well. 
This important problem should be resolved experimen- 
tally in the nearest future. Finally, further violation of 



the KNO-scaling in multiplicity distributions is demon- 
strated in QGSM. The origin of both conservation and 
violation of the scaling trends is traced to short range 
correlations of particles in the strings and interplay be- 
tween the multi-Pomeron processes at ultra-relativistic 
energies. 
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